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$_i \ Abstract 
Q+ 

The next to leading order (NLO) contribution of the generalized finite size mass 
£f) ' shift formula for an interacting two stable particle system in a periodic L 3 box is 

discriminated with maintaining its model independent structure and validity to all 
orders in perturbation theory. The influence of the NLO contribution is examined 
' for the nucleon mass shift in the realistic nucleon-pion system. 
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rin [ Measurements of hadron spectrum in unquenched lattice QCD simulations al- 

ways suffer the finite volume effect from the associated virtual cloud of lightest 
particles in the spectra, which may lap the whole lattice once or several times 
owing to periodic boundary conditions. Finite size mass shift formulae, in- 
volving the quantum loop effect of pions in finite volume, are thus derived 
in chiral perturbation theory (ChPT) [1,2,3,4,5,6,7,8] and applied to the sim- 
ulation results for the purpose of controlling its volume dependence and of 
identifying the value corresponding to the thermodynamic limit (see [9] for a 
recent review). 

In our previous paper, we looked at this issue [10] from a general field theoret- 
ical point of view (without sticking to ChPT) and derived the finite size mass 
shift formula for the interacting two stable particle system in a periodic L 3 
box, as an extension of Liischer's formula for self-interacting bosons [11,12]. 
Remarkable points of Liischer's formula are that the finite size mass shift in 
a periodic box is related to forward elastic scattering amplitudes in infinite 
volume, which is model independent, and can be valid to all orders in pertur- 
bation theory up to a certain error term [12]. 

In perturbation theory the physical mass is denned from the pole position of 
the full propagator. Using this fact, the finite size mass shift of a (bosonic) 
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particle 1 Ara(L) = M(L) — rain Euclidean space can be defined as 

Am(L) = -^-[S L (p)-S(p)]+0((Am) 2 ) at p=(im,6), (1) 

where an d denote the self energies in the finite and infinite 

volumes, respectively. We renormalize the self energy in infinite volume as 
= a^2^G°) = at p 2 = — m 2 . In perturbation theory contains the 

number of sums over discrete spatial loop momenta, cfi(L) = 2ixrti/ L (n^ e Z 3 ), 
depending on the number of loops (i = 1, . . . , N\ oop ). These summations can 
be rewritten as integrals by using the Poisson summation formula. Then the 
integrand of is reduced to the same form as that of E(p) apart from the 

exponential factors e~ tLmi '^ and summations over integer vectors m, G Z 3 . 

Since the magnitude of rhi e Z 3 acts as the weight of the exponential suppres- 
sion factor of the mass shift formula, the leading order (LO) contribution to 
Am(L) for an asymptotically large L can be specified by requiring that only 
one of them has a non-zero value |m| = \rhi\ = 1 and the others have zero 
\rhj\ = (j i). In other words, the asymptotic formula can be described 
by the collection of the effective one-loop diagrams with an exponential factor 
e -iLm-q^ w h ere the other loop integrals without exponential factors are reduced 
to the part of the definition of the vertex function in infinite volume. Luscher 
originally discussed this case [11,12] and we also did it in the previous pa- 
per [10]. The order of the error term in the formula was then consistent with 
that of the next to leading order (NLO) contribution; \rh\ = \rhi\ = \[2 and 

For the realistic application of the formula to analyzing lattice data, however, 
it is desirable to reduce the ambiguity associated with the error term. In the 
present paper, we thus aim to discriminate the NLO contribution (|m| = \/2) 
in the formula for the two particle (A-B) system, in particular, while maintain- 
ing its model independent structure and validity to all orders in perturbation 
theory. In this case, the task is still the same as for the \rh\ — 1 case; we eval- 
uate the effective one-loop diagrams as listed in Fig. 1. We may here assume 
that A particle carries a conserved charge, so that interaction induced by the 
three-point vertex AAB and four-point charge conserving vertices are taken 
into account. It should be noted that what is nontrivial for such an extension 
is not so much evaluating the \rh\ = \/2 contribution itself as evaluating the 
\m\ = 1 contribution with an error term at most of the order of the NNLO 
contribution (|m| = \/3). Otherwise the NLO contribution will be obscured 
in the error term. In fact, it is straightforward to compute the \rh\ = \[2 
contribution once the procedure is established for |m| = 1. 

1 The fermionic mass shift can also be defined in a similar way by sandwiching the 
self energies between spinors. 
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Fig. 1. Effective one- loop self-energy diagrams which contribute to the mass shift 
formula in the bosonic A-B system. Solid lines with an empty circle correspond 
to the propagator of A particle, Ga, and dashed lines to that of B particle, G B - 
Shaded blobs are vertex functions, Taab, ^aaaa, ^aabb, ^bbb, at certain orders 
in perturbation theory. It is assumed that A particle carries a conserved charge. 

The result turns out that for the mass ratio 
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with a max w 



0.418, it is possible to discriminate the NLO contribution and 
the final expression can be written as 
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where the first and second lines correspond to the LO and NLO contributions, 
respectively. In this expression, 
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and F AB (v) denotes the forward scattering amplitude of A(p) + B(q) — > A(p) + 
B(q) in infinite volume (is = iqo the crossing variable). F AB (is) has poles at 
is = ±is B . The coupling A is then defined by exploiting the residue of F AB (is) 
as 



\i^y-vl)F AB {v) = ^. (5) 



The basic line of the derivation of Eq. (3) is the same as in our previous 
paper [10], where the detailed notation of the propagators and the vertex 
functions are also given. In what follows, we shall present a derivation for the 
case that both A and B particles are bosons. The extension to the case that 
particle A is a fermion is straightforward, and the final result is exactly the 
same as in Eq. (3), which is one of the advantages of the model independence 
of the formula. Here, we concentrate on evaluating the diagram (bl) in Fig. 1, 
which is typical for the two particle system and provides us with a key idea 
how to control the error term. The other diagrams are then evaluated in a 
similar way. 

The self-energy diagram (bl) for \m\ = 1 is expressed as 



4f l=1) =6 / -^-e iL ^r AAB (-p,(l-r ] )p + q ;r ] p- q )G A ((l-r ] )p + q ) 



xG B (rjp - q)T AAB (p, -(1 - rf)p - q; -rjp + q) 



_ , (6) 

p={im,A, 0) 



where 77 is a real parameter at least in the range [0, 25(a)] for 5(a) e (0, 1/2]. 
For the purposes of evaluating the integral i] can be chosen appropriately 
depending on the mass ratio a. Our concern is whether there exists such a set 
of i] and 5 for a given a to make the error terms smaller than the desired order 
of magnitude. For our purpose this is 0(e~ L7M ) with 7 = \/3. In our previous 
work [10], we chose r] — 5 — a 2 /2, which was sufficient to control the error 
term up to 0(e _v ^ iM ). But this choice turns out to be inappropriate in the 
present case (see our final choice in Eq. (12)). The overall factor 6 originates 
from rotational invariance among q±, q-2 and ^3. 

Firstly, we perform the complex q x contour integration by focusing on the 
poles of the propagators G A and G B of one-particle states at 2 



We are assuming that there is no bound state below the two-particle threshold. 
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q[ A) = i^ql + q\ + (2 V - r?)m\ + i2(l - r ? )m A g , (7) 
q[ B) = i\J ql + q\ + m 2 B - r] 2 m 2 A - i2r)m A qo , (8) 

in the complex qi upper half plane, respectively, 3 where q\ = q\ + q\- We 
may set a contour which goes along the real q\ line and the line Im q\ = Q\ > 
closed at ±oo to pick up the residues at q[ A ^ and/or q[ B \ In order to relate 
the mass shift to the on-shell forward scattering amplitude like in Eq. (3), we 
find at this point that the upper path 6\ must be chosen so as to satisfy the 
following conditions; 

(i) the contribution from the upper path itself is smaller than the error term 
0(e" Lw ), 

(ii) the contour covers the range of Im q[ A ^ and/or Im q[ B ^ for q and q± in a 
certain ball 

M = {(q ,q ± )ER 3 \q 2 + ql<v 2 } , (10) 

(iii) the contour picks up no residue except for the poles at q[ A ^ and/or q[ B \ 

Here the condition (iii) must be guaranteed even if q is extended to a complex 
variable and shifted as qo ^ Qo — «(1 _ i]) m A for q[ A ^ and/or q — > q + irjuiA 
for q[ B \ where q satisfies the on-shell condition q 2 = —m 2 A and/or q 2 = —m 2 B . 

To examine the condition (iii), we use the fact that the vertex function 
Faab(TP, (1 - v)p ^ Q'j VP -F q) with r\ e [0, 25], initially defined for (p, q) e 
K 4 x R 4 , is analytic inside the complex domain 



D = {(p,?)GC 4 xC 4 (Tm{(l- 5) p±^q}^j < m\ , 

(lm{5p±^q}^ <m|} . (11) 

The basic observation for finding this domain is that the vertex function at 
any higher order in perturbation theory consists of a set of A and B lines (free 
propagators) . The denominator of the Zth A or B line is then parametrized as 
(k(l) +r(/)) 2 +m 2 A or (k(l) +r(l)) 2 + m 2 B , where k(l) is the external momentum 
flow given by a combination of complex variables p and q, and r(l) is a combi- 
nation of internal loop momenta to be integrated out, which is a real variable 
in Euclidean space. It then follows that the vertex function has no singularity 

3 For a,kR, 

Im (iV^Tib) = \J(Va 2 + b 2 + a)/2 . (9) 
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Fig. 2. The external momentum flow in the AAB vertex function. Arrows represent 
the flow direction. 

if (Im k(l)) 2 < m\ and (Im k(l)) 2 < m 2 B for all A and B lines. In order to find 
the possible choices of k(l), we may label the three bare vertices where the 
external momenta, p, (1 — rf)p + q and r/p — q, are plugged in (and out) as a±, 
0,2 and b, respectively (e.g. a\ = = b at the tree level). Note that whenever 
A particle carries a conserved charge, there always exists a set of A lines con- 
necting di and a,2- In Fig. 2, we show the possible external momentum flow 
inside the AAB vertex function; they are basically classified into two cases, 
the connected A lines flow through b (left) and the connected A lines do not 
flow through b (right). One can add any internal lines depending on the order 
in perturbation theory, which however carry no external momentum and do 
not affect the singularity of the vertex function. Inserting the largest exter- 
nal momenta for A and B lines into (Im k(l)) 2 < m\ and (Im k(l)) 2 < m 2 B , 
respectively, one can specify the domain O as in (11). 

We then realize that there is no integration path 0\ at any value of a which 
satisfies all conditions (i)~(iii) for both poles q[ A ^ and q± simultaneously. 
However, we find that it is possible to choose 6\ = ^ "f 2 /j, 2 + rf"m 2 A with the 
upper bound of the ball v 2 = 7 2 /i 2 — m 2 B + r] 2 m 2 A in Eq (10), which satisfies 

( B) 

the conditions only for q[ within a limited range of a. The choice of rj and 
5 is quite subtle, but by choosing 



r] — c 8 , 5 = a 

with c = 0.95, the allowed range of a is maximized as a e (0, a max ], where 
ttmax ~ 0.418. 4 Thus we obtain 



4 In this region the contribution from q[ A) can be neglected since it is of 0(e 

To show this explicitly, we need to carry out the qo contour integration by performing 

(A) 

the momentum shift q$ — > qo + iQ^rriA for q\ . Here, the domain D constrains 9q 
to be 6 < ((4 - c)S - 2(1 - c 2 /A)5 2 )/{l - (2 - c)<5), 9 > -((4 - c)S - 2(1 - 
c 2 /4)5 2 )/(3 - (2 + c)<5), > -(2a 2 - cS - 2(1 - c 2 /4)5 2 )/(l + (2 - c)<5), and 9 > 



4~7 2 (l-f) 
2(c 2 + 2c + 2) 



(12) 
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Fig. 3. qo integration contour. 




6i / fo \s ^ — r ^B(-p, (i - v)p + q;vp- q) 

Jb (Ztt) 4 Zqi 
xG A ((l - r])p + q)T AAB (p, -(1 - r])p - q; -rjp + q)\ 



«i=«i 



(13) 



We remark that if we set 7 = 2, where 0(e~ 2L ^) corresponds to the NNNLO 
contribution, there is no parameter set rj, 5, v) which fulfills the above 
conditions at any value of a, indicating that one cannot discriminate the 
NNLO contributions along this line. 

Secondly, we perform the complex q contour integration along the path in 
Fig. 3. The path CD is parametrized by shifting the momentum q — > q + 
irjrriA- On this path the argument of G A becomes p + q, where q satisfies the 



on-shell condition q 2 = —m 2 B , since q[ = i\J ql + q\ + rn 2 B . Inside the contour 

the integrand has no pole except at 



which is guaranteed by the above condition (iii). Note that < Im q' < rjm A 
for the value of rj specified by Eq. (12). The contributions from the paths BC 
and DA are at most 0(e _L7M ) since Im q[ B ^ > 7//, which is due to the choice 
of the ball B with v 2 = ^ 2 /j, 2 — m 2 B + rfm 2 A in (10). Thus the integral (13) is 
replaced by one along the path CD with the residue contribution at q' 



- (2a 2 -cS- 2(1 - c 2 /A)5 2 ) / (l - (2 + c )S). On the other hand, we find Im q\ ' > jfi, 
if 9q < 7] — 1 + y/l — 7 2 a 2 (l — a 2 /4). Solving these inequalities (numerically) with 
Eq. (12), we find that c = 0.95 yields the maximal value of a. 





(14) 
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Fig. 4. Ingredients of F A b(v) in the system. The labels represent the corre- 
spondence with the self-energy diagrams in Fig. 1. 



(|m|=l)_o 

1 bi — 6 



d 2 q± e-^v 7 ^ A 2 



' (2tt) 2 2y/<g_ + y 2 2 "b 

+ 6 f ^ V ^ ql+ql+m% F <*-V (iQn ) + 0(e-^) (15) 
+ 6 i (2*)" 2 ^ + gi + m | F - ( ^ 0) + ° (e } ' (15) 



where 



f ab 1 \i%)= r AAB(-p,P + q; -q)G A (p + q) 

xT A ab(p, -v - q; q) \^=-m\, q 2 =-m% ( 16 ) 

corresponds to a one-particle-irreducible (1PI) part of the forward scattering 
amplitudes of F AB (y = iq ), graphically represented as (bl-1) in Fig. 4. By 
using the crossing relation Ff^f^^v) = F A b B ~ 2 \v), one can replace Ff^f^^u) 
in Eq. (15) by {F A V ^~ 1 \v) + F A b B ~ 2) (v)) /2. In the first term in Eq. (15), the 
effective renormalized coupling A is defined by Eq. (5), where F AB ^{y) or 
F AB ~ 2 \u) has the pole at v — ±u B , and the integral region is specified by 
inserting q = q' to B: 



a 2 

V = {q ± G M 2 | q 2 ± < ( 7 2 - % 2 + 2 V m 2 A (v - y)} ■ (17) 



We then carry out the q± integration in Eq. (15) by using the integral formula 
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-L P 



-oo (2vr) 2 2v /g2 + p 2 AttL 



(18) 



where the integral region can be extended from B or B' to infinity, because 
the boundary contributions of B and B' are already smaller than the order of 
the error term. Hence, we end up with 



A\m\=l) 
'61 



1 /6 

8^ \L 



2u B 
+0{e~^) . 



(19) 



Other self-energy diagrams in Fig. 1 can be evaluated in a similar way up 
to 0(e~^ L ^), yielding the corresponding 1PI part of the forward scattering 
amplitude. Note that the contributions from the self-energy diagrams (al) and 
(cl) are already smaller than 0(e~^ Lfl ) for a e (0,a max ]. By combining all 
contributions we can discriminate the \rh\ = 1 contribution up to 0(e _v ^ L/i ). 

The |m| = \f 7 2 contribution is given by the integral 



/f M) = 12 / ^ 4 e iL ^r AAB (-p, (1 - V )p + q; VP - q)G A ((l - V )p + q) 



(2tt)' 

xG B (r}p- q)T AAB (p, 



[l-r])p-q; -rjp + q) 



(20) 



p=(im A , 0) 



Rotating the q\-q 2 axis by 7r/4, we define q± = (qi + q2)j\[2 and q 2 = (—qi + 
q2)j\p2. Then, apart from the new exponential factor e 1 ^^ 1 and an overall 
factor 12, the integrand becomes exactly the same as in Eq. (6). Thus the 
evaluation is straightforward and the result is 



A\m\=V2) _ J_ / 12 \ 



2u B 



+ r d ^e-^V^ {F ^\iq ) + F 

J-ogA'K 



(»-2)/.- 
AB 



{no)} 



+o( e -^) , 



(21) 



where the error term becomes automatically smaller than the \m\ = 1 case. 
Evaluating other diagrams similarly and combining the result of |m| = 1, we 
arrive at the mass shift formula in Eq. (3). 



Finally, let us examine the influence of the \m\ = \/2 contribution by looking 
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0.1 




Fig. 5. The nucleoli mass shift as a function of £ = Lm^. 

at the nucleon mass shift in the realistic iV-7r system with = 938 MeV 
and m n = 140 MeV. As the mass ratio is a = m^/m^ = 0.149, Eq. (3) is 
applicable. Moreover, since the formula is valid to all orders in perturbation 
theory and is expected to hold nonperturbatively, we may insert the empirical 
N-Ti scattering amplitude into Eq. (3). The subthreshold expansion of the N-tt 
forward scattering amplitude around v = is parametrized as [13] 

F Nv (v) = 6rn N D+{u) , (22) 
where 



D » = ^TTX^ + d ™ m ^ + m - V + 4) + 0{v & ) • (23) 

Til U g £/ 

The isospin sum is taken in Eq. (22), neglecting the effect of isospin symmetry 
breaking. The coupling constant is g 2 = 14.3. The first term in Eq. (23) 
is identified with the pseudovector nucleon Born term with z/g = m^/2mjv m 
0.07m n . The effective coupling is then easily computed by using Eq. (5) as A 2 = 
—Ylg 1 v\ 3 . The coefficients of the other terms are given by dgo = —1-46(10), 
df = 1.12(2) and d% = 0.200(5) [13]. Hereafter we only take into account the 
mean of these values. 

In Fig. 5, we plot A(£ = Lm^) = Am^/inN for the |m| = 1 and |m| = y/2 
contributions (dotted and dashed lines, respectively) and the sum of these 
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contributions as \m\ < a/2 (solid line), where £ = 1 corresponds to L = 1.4 fm. 
It reveals that the \m\ = a/2 contribution is quite large for the plotted region 
of £. For instance at £ = 2 the mass shift is expected to occur more than 
1.2 (\m\ = 1) + 0.8 (\m\ = a/2) = 2.0 % (> 20 MeV). If one estimates the 
|m| = a/3 contribution itself at £ = 2, 5 this merely contributes the mass shift 
by 0.2 %. This is due to the smaller geometrical factor (e.g. 6 for \rh\ = 1, 
12/ V2 « 8.49 for |m| = a/2, and 8/a/3 » 4.62 for |m| = a/3) as well as the 
larger exponential decay factor. In this sense the nucleon mass shift formula 
is significantly modified by discriminating the NLO contribution. Note that 
the negative mass shift for £ < 1 is due to the contribution from the term 
involving the N-ir forward scattering amplitude (ingredients of the |m| = 1 
curve can be found in Ref. [10]). 

To summarize, we have investigated the finite size mass shift formula for the 
two stable particle system in a periodic L 3 box. We have found that it is pos- 
sible for the mass ratio a G (0, a max ] with a m ax ~ 0.418 to discriminate the 
NLO contribution with maintaining its model independent structure and va- 
lidity to all orders in perturbation theory. The final expression is then written 
as in Eq. (3). Along the way we have also realized that it is impossible to 
discriminate the NNLO contribution along the line discussed above once the 
error term is set by 7 = 2. In fact, in order to discriminate more higher order 
contributions, one should go back to the definition of the mass shift in Eq. (1). 

We are grateful to the members of lattice forum in DESY theory group in 
Hamburg, in particular, H. Wittig and I. Montvay for valuable discussions 
and comments. We appreciate fruitful comments from P. Weisz. Y.K. thanks 
T.R. Hemmert for useful discussions at the meeting of the DFG Forscher- 
gruppe 'Lattice Hadron Phenomenology' at DESY-Zeuthen in February. 
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